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I. INTRODUCTION 
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In this work we consider the 1-component real scalar 
(f> 4 theory in 4 space-time dimensions on the lattice and 
investigate the finite size scaling of thermodynamic quan- 
tities to study whether the thermodynamic limit is at- 
tained. The results are obtained for the symmetric phase 
of the theory. 



II. THERMODYNAMIC SET UP 



The thermal expectation value of an operator O is 



(1) 



where Z((3) is the partition function. For our finite tem- 
perature study, we have kept the temporal side fixed 
at T — /3 _1 = (7V t a t ) _1 and varied the spatial extent 
L = N s a s by going to higher values of LT to exam- 
ine the thermodynamic limit. Energy density E and 
pressure P are obtained by the following derivatives 
of the partition function: E = T 2 /V d\nZ/dT\ v = 
-l/{N^a 3 s N t ) d\nZ/da t \ v and P = T d\nZ/dV\ T = 
l/(3N 3 alN t a t ) d\nZ /da s \ T where a s and at are the lat- 
tice spacings along the spatial and temporal directions. 
The continuum Euclidean action for the scalar </> 4 theory 



S[<(>] = J dtj d\^d^(x)d^(x) + \ml<t>{xf + ^<f>(x) 

(2) 

The following redefinitions are made to trade off the 
set of couplings (go, mo) for (k, A): 
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(3) 



where £ = a s /at is the anisotropy parameter. This yields 
the conventional form of the action for doing simulations: 
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S[4>] = y~l{~ K s 4>{xa + ai)<f>(x a ) - K t (j)(x a + a t )(j>(x a ) 

X a 1 — 1 

+4>{x a f + \(4>{x a ) 2 ~ l) 2 } (4) 

Since our simulations are on lattices with £ = 1, 
we calculate the relevant thermodynamic expressions for 
anisotropic lattices and then set the spacings to be equal. 
We renormalise the pressure to be zero at T = by sub- 
tracting the corresponding zero temperature operators 
from their finite temperature counterparts. This proce- 
dure is adopted to renormalise all the operators appear- 
ing in the other thermodynamic quantities as well and 
is denoted by adding the subscript "subt" to the opera- 
tors. Finally, the expressions for the energy density E, 
pressure P and E — 3P in this formulation is: 
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(5) 



where bracketed quantities denote thermal averages. 
Note that the quantities Ea\, Pa\ and (E—3P)at are all 
dimensionless quantities and will henceforth be referred 
to as E, P and E — 3P respectively in lattice units. The 
derivatives of the couplings of the theory with respect to 
£ are the anisotropy coefficients. 

The renormalisation conditions for the theory involve fix- 
ing the 2- and 4-point renormalised vertex functions at 
zero momenta as follows: 



3t = -rg>(0); «te = -1^(0,0,0,0) 



( 4 )f 



(6) 



Equation ([5]) requires evaluation of the following 
derivatives: 



«k s = at q — 
oa t 



, dn t 

O-Kt = a t-^— 

oa t 



; ax = at — 
^ oa t 



(7) 



cx = A'| 0t (8) 



where /' = df/d£. The first three derivatives can 
be completely specified using [3 = at dgo/dat\^ and 

7 = at drfiQ 2 /9a* L functions. We have used the 1- 
loop expressions for the j3 and 7 functions of the lat- 
tice theory available in literature [![. The latter three 
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FIG. 1: The Feynman diagrams to 1-loop order that need 
to be calculated. The first set of diagrams is for the 2-point 
vertex function while the second one is for the 4-point vertex 
function. 



involve the anisotropy coefficients g' Q = dgo/d£,\ at and 
(too 2 )' = <9too 2 /<9£| , whose expressions do not exist in 

I (If 

literature for this scalar theory. We have evaluated the 
anisotropy coefficients upto 1 loop order by summing the 
Feynman diagrams shown in fig [1] to be: 



K{x) 
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[(3m| + 12)/ 4 (2a;) - 12/g(2a:)Ii(2a:)] 
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x 2 dxe- x{ ™ 2 « +s) K{x) 



(9) 



where /„ (x) denotes the modified Bessel function of order 
n. The calculation of these coefficients was done using 
Mathematica which could handle large range of numbers, 
much beyond the capacity of usual C compilers. This was 
specially important since the numerical values depend 
on a delicate cancellation among Bessel functions. The 
anisotropy coefficients evaluated at the simulation points 
is listed in table |TJ 



III. NUMERICAL METHODS 

A overrelaxation+Metropolis (OR+M) algorithm was 
used to generate the equilibrium configurations on which 
the operators were measured. The OR+M algorithm 
caused the correlation time at critical point to reduce 
by a factor of 10 compared to the Metropolis and the 
time to generate independent configurations was 0.4s. 
Configurations separated by 2>T int were considered to be 
independent, where Tint was the Monte-Carlo time for 
the autocorrelation function to decay to 1/e. To check 
the accuracy of our program, simulations were done at 
(A = 0.275376, k = 0.288998) where we measured the 
vacuum expectation value of the 4> field and its suscepti- 
bility, and compared it with the results obtained in [2|. 

The jackknife method was used to calculate the errors 
on the operators. It was checked that the block size used 
was much above the minimum bin size over which the 
variance showed a plateau. 



TABLE I: Parameter values at each of the different points, gn 
is calculated in 1 loop perturbation theory. For compari- 
son, estimate of mass in 1 loop perturbation theory is also 
shown. 



JNo. 


A 


T ) 

D 




1 1 

L) 


K 


0.2516 


0.2526 


0.253452 


0.25438 


A 


0.011153 


0.0117 


0.0122 


0.0134 


9r 


3.1772 


3.18955 


3.16478 


3.17108 


ma 


0.2853 
±0.0002 


0.2376 
±0.0003 


0.1897 
±0.0002 


0.1416 
±0.0002 


ma (in 1 loop) 


0.3012 


0.2586 


0.2172 


0.1834 


g'o 


-1.9372 


-2.06348 


-2.16451 


-2.34449 


(mo 2 )' 


0.367901 


0.369738 


0.367193 


0.368173 



Subtraction dependence in point A runs at T/m=0.877 
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10 15 20 25 30 35 40 
Lattice used for T=0 subtraction 



8°X4 * 12°X4 ; □ : 16 X4 • 

Subtraction dependence in point A runs at T/m=0.877 
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Subtraction dependence in point A runs at T/m=0.877 
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FIG. 2: Comparison of different subtraction schemes on the 
E/Esb, P/Psb and (E - 3P)/T 4 at point A 
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Comparison between the N t =4 and N,=6 results 
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Comparison between the N t =4 and N,=6 results 
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Comparison between the N t =4 and N,=6 resuits 
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FIG. 3: Scaling of the E/T 4 , P/T 4 and (E - 3P)/T 4 for 
Nt = 4 and N t = 6. 



IV. PHYSICS RESULTS 

The simulations were done on Renormalisation Group 
(RG) trajectories of constant gR. A point was chosen 
(k = 0.2516, A = 0.011153) where the scaling violations 
could be neglected [3j. The physical mass was extracted 
non-perturbatively at this point A. Points B, C and D 
with cut-off effects 4/5, 2/3 and 1/2 of point A and lying 
on the same RG trajectory were first constructed within 
1-loop perturbation theory and then non-perturbatively 
tuned to get the desired ma. The parameter details of 
the points are given in table HI 



To estimate zero T subtraction, operators evaluated on 
lattices ranging from 40 4 till Nf were subtracted from the 
corresponding operators evaluated on the x 4 lattice 
at point A. The result is shown in fig [21 Our results 
indicate that a plateau is reached for all the thermody- 
namic quantities when the subtraction lattices used iV 4 
are much greater than TV 4 except for the runs at 32 4 and 
36 4 which seem to give erratic results. We checked that 
the problem was not in the analysis. In what follows, 
we do not include these two sets in the discussion of the 
thermodynamic limit. The result of subtracting with a 
40 4 lattice agrees with the results obtained by subtrac- 
tion from the N 4 " lattices except when N is comparable 
to N s . Thus, we chose to subtract using the largest lat- 
tices that we have simulated which is 40 4 . 
The E, P and E - 3P are plotted vs LT in fig GO To 
establish the thermodynamic limit at constant physics 
we compare the LT dependence of the Nt — 4 and the 
N — 6 lattices at each of the four different points. Our 
results show that as we decrease ma at fixed LT the finite 
size effects become larger, as expected. At fixed m and 
for the same N t , the absolute value of the slope decreases 
with larger LT. In particular, for the N t — 6 lattices, for 
energy density and pressure, the local slope between the 
LT = 3 and LT = 4 is zero within error. The Nt — 4 and 
N = 6 lattices at fixed m can be compared by looking at 
the slopes of the thermodynamic quantities as a function 
of LT. The local slope for N t = 6 would be lesser than 
the Nt = 4 counterparts for the same range of LT. This 
comparison also yields the expected results except in the 
case of the smallest mass. 

For E and P, the thermodynamic limit is reached for 
the LT — 4 runs for the masses ma — 0.2853 and 
ma = 0.2376, since the local slope is zero within error. 
Even though this is not the case for E — 3P, it certainly 
is close to the thermodynamic limit since the values ap- 
proach a plateau. For example, at point A, the E — 3P 
for Nt — 6 lattice decreases from (6.16 ± 0.02) x 10~ D 
at LT = 2 to (4.212 ± 0.006) x 10~ 5 at LT — 3 and to 
(4.025 ± 0.004) x 10" 5 at LT = 4. For the last two 
points where the finite size effects are large the situa- 
tion is less clear. In terms of ml, our results show that 
mL = 4.55 is not enough for the thermodynamic limit, 
while by mL — 5.70 the limit is attained. We conclude 
that to reach the thermodynamic limit, one needs LT ~ 4 
and mL > 4.5. We also found that the T = subtraction 
should be made with lattices with N sufficiently greater 
than N s . In this work we have not tried to take the 
continuum limit. 
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